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Abstract. We establish the existence, stability, and asymptotic behavior of 
transonic flows with a transonic shock past a curved wedge for the steady 
full Euler equations in an important physical regime, which form a nonlinear 
system of mixed-composite hyperbolic-elliptic type. To achieve this, we first 
employ the transformation from Eulerian to Lagrangian coordinates and then 
exploit one of the new equations to identify a potential function in Lagrangian 
coordinates. By capturing the conservation properties of the system, we derive 
a single second-order nonlinear elliptic equation for the potential function in 
the subsonic region so that the transonic shock problem is reformulated as a 
one-phase free boundary problem for the nonlinear equation with the shock- 
front as a free boundary. One of the advantages of this approach is that, given 
the shock location or equivalently the entropy function along the shock-front 
downstream, all the physical variables can be expressed as functions of the 
gradient of the potential function, and the downstream asymptotic behavior of 
the potential function at infinity can be uniquely determined with a uniform 
decay rate. To solve the free boundary problem, we employ the hodograph 
transformation to transfer the free boundary to a fixed boundary, while keeping 
the ellipticity of the nonlinear equation, and then update the entropy function 
to prove that the updating map has a fixed point. Another advantage in our 
analysis is in the context of the full Euler equations so that the Bernoulli 
constant is allowed to change for different fluid trajectories. 
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1. Introduction. We are concerned with the existence, stability, and asymptotic 
behavior of steady transonic flows with transonic shocks past curved wedges for the 
full Euler equations. The two-dimensional steady, full Euler equations for polytropic 
gases have the form (c/. [14, 15, 25]): 

( V • (pu) = 0, 

<1 V • (pu0u)-t Vp = 0, (1.1) 

[ v-(pu(£; + H))=o, 

where V = Vx is the gradient in x = (xi, X 2 ) € u = (ui, U 2 ) the velocity, p the 
density, p the pressure, and E = Ijup -|- e the total energy with internal energy e. 

Choose density p and entropy S as the independent thermodynamical variables. 
Then the constitutive relations can be written as (e,p, 9) = (e(p, S),p{p, S),9{p, S)) 
governed by 

9dS = de — ^dp, 

P 

where 9 represents the temperature. Eor a polytropic gas, 

R 

p = (7 - l)pe, e = Cv9, 7 = 1-1-> 1, (1.2) 

Cy 

or equivalently, 

p = p{p,S) = e = 6(^,5) =(1.3) 

7-1 

where i? > 0 may be taken to be the universal gas constant divided by the effective 
molecular weight of the particular gas, c„ > 0 is the specific heat at constant volume, 
7 > 1 is the adiabatic exponent, and k > 0 is any constant under scaling. 

The sonic speed of the flow for polytropic gas is 



The flow is subsonic if juj < c and supersonic if juj > c. Eor a transonic flow, both 
cases occur in the flow, and then system (1.1) is of mixed-composite hyperbolic- 
elliptic type, which consists of two equations of mixed elliptic-hyperbolic type and 
two equations of hyperbolic type. 

System (1.1) is a prototype of general nonlinear systems of conservation laws: 

Vx-r(f7)=o, xeR”, (1.5) 

where U : R" R™ is unknown, while F : R™ —>■ is a given nonlinear 

mapping for the mxn matrix space For (1.1), we may choose U = (u,p, p). 

The systems with form (1.5) often govern time-independent solutions for multidi¬ 
mensional quasilinear hyperbolic systems of conservation laws; cf. Lax [20, 21]. 

It is well known that, for a steady upstream uniform supersonic flow past a 
straight-sided wedge whose vertex angle is less than the critical angle, there ex¬ 
ists a shock-front emanating from the wedge vertex so that the downstream state 
is either subsonic or supersonic, depending on the downstream asymptotic condi¬ 
tion at infinity (see Appendix B and Fig. 2.1 for the shock polar). The study of 
two-dimensional steady uniform supersonic flows past a straight-sided wedge can 
date back to the 1940s {cf. Courant-Friedrichs [15]; also see Prandtl [26] and von 
Neumann [28]). 

For the case of supersonic-supersonic shock {i.e. both states of the shock are 
supersonic), local solutions around the curved wedge vertex were first constructed 
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by Gu [19], Li [23], Schaeffer [27], and the references cited therein. Global potential 
solutions are constructed in [6, 7, 15, 29] when the wedge has certain convexity or 
the wedge is a small perturbation of the straight-sided wedge with fast decay in the 
flow direction. Furthermore, the stability and uniqueness of entropy solutions in 
BV containing the strong supersonic-supersonic shock were established in Ghen-Li 

[4]. 

For the case of supersonic-subsonic shock (z.e. transonic shock-front), the sta¬ 
bility of these fronts under a perturbation of the upstream flow, or a perturbation 
of wedge boundary, has been studied in Chen-Fang [11] for the potential flow and 
in Fang [16] for the Euler flow with a uniform Bernoulli constant. In particular, 
the stability of transonic shocks in the steady Euler flows with a uniform Bernoulli 
constant was first established in the weighted Sobolev norms in Fang [16], even 
though the downstream asymptotic decay rate of the shock slope at infinity was 
not derived. 

In this paper, one of our main objectives is to deal with the asymptotic behavior 
of steady transonic flows with a transonic shock past a curved wedge for the full 
Euler equations, especially the uniform decay rate of the transonic shock slope and 
the subsonic flows downstream at infinity. For a fixed uniform supersonic state , 
there is an arc on the shock polar corresponding to the subsonic states; see Fig. 
2.1. When the wedge angle is less than the critical angle 0^, the tangential point T 
corresponding to the critical angle divides arc HS into the two open arcs TS and 
TH. The nature of these two cases is very different. In this paper, we focus mainly 
on the stability of transonic shocks in the important physical regime TS when the 
wedge angle is between the sonic angle 0^ and the critical angle 0^ > 0w- 

To achieve this, we first rewrite the problem in Lagrangian coordinates so that the 
original streamlines in Eulerian coordinates become straight lines and the curved 
wedge boundary in Eulerian coordinates becomes a horizontal ray in Lagrangian 
coordinates. Then we exploit one of the new equations to identify a potential 
function (j) in Lagrangian coordinates. By capturing the conservation properties of 
the Euler system, we derive a single second-order nonlinear elliptic equation for the 
potential function </> in the subsonic region as in [1], so that the original transonic 
shock problem is reformulated as a one-phase free boundary problem for a second- 
order nonlinear elliptic equation with the shock-front as a free boundary. One of 
the advantages of this approach is that, given the location of the shock-front, or 
equivalently the entropy function A (which is constant along the fluid trajectories) 
along the shock-front downstream, all the physical variables U = {u,p,p) can be 
expressed as functions of the gradient of (j), and the asymptotic behavior (j)°° of the 
potential (j) at the infinite exit can be uniquely determined. 

To solve the free boundary problem, we have to determine the free boundary, and 
both the subsonic phase and entropy function defined in the downstream domain 
with the free boundary as a part of its boundary. We approach this problem by 
employing the hodograph transformation to transfer the free boundary to a fixed 
boundary, while keeping the ellipticity of the second order partial differential equa¬ 
tions, and then by updating the entropy function A to prove that the updating map 
for A has a fixed point. 

For given entropy function A, we first determine a priori the limit function of 
the potential function downstream at infinity. Then we solve the second-order 
elliptic equations for the potential function in the unbounded domain with the fixed 
boundary conditions and the downstream asymptotic condition at infinity. This is 
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achieved through the fixed point argument by designing an appropriate map. In 
order to define this map, we first linearize the second-order elliptic equation for 
the identified potential function </> based on the limit function (j)°° of (j), solve the 
linearized problem in the fixed region, and then make delicate estimates of the 
solutions, especially the corner singularity near the intersection between the fixed 
shock-front and the wedge boundary. These estimates allow us to prove that the 
map has a fixed point that is the subsonic solution in the downstream domain. 
Finally, we prove that the entropy function A is a fixed point via the implicit 
function theorem. 

Since the transformation between the Eulerian and Lagrangian coordinates is in¬ 
vertible, we obtain the existence and uniqueness of solutions of the wedge problem in 
Eulerian coordinates by transforming back the solutions in Lagrangian coordinates, 
which are the real subsonic phase for the free boundary problem. The asymptotic 
behavior of solutions at the infinite exit is also clarified. The stability of transonic 
shocks and corresponding transonic flows is also established by both employing the 
transformation from Eulerian to Lagrangian coordinates and developing careful, 
detailed estimates of the solutions. 

Another advantage in our analysis here is in the context of the real full Eu¬ 
ler equations so that the solutions do not necessarily obey Bernoulli’s law with a 
uniform Bernoulli constant, i.e., the Bernoulli constant is allowed to change for 
different fluid trajectories (in comparison with the setup in [9, 10, 13, 16]). 

By the closeness assumption of solution U to the uniform flow in the subsonic 
region, we obtain the asymptotic behavior of 1/ as j/i —>■ oo. The asymptotic state 
U°° = {u°° ,p°°, p°°) is uniquely determined by state U~ of the incoming flow and 
the wedge angle at infinity. 

We remark that, when Uq is on arc TH (see Fig. 2.1 below), the nature of 
the oblique boundary condition near the origin is significantly different from the 
case when Uq is on arc TS. Such a difference affects the regularity of solutions 
at the origin in general. It requires a further understanding of global features of 
the problem, especially the global relation between the regularity near the origin 
and the decay of solutions at infinity, to ensure the existence of a (7^’“ solution. A 
different approach may be required to handle this case, which is currently under 
investigation. 

The organization of this paper is as follows: In §2, we first formulate the wedge 
problem into a free boundary problem and state the main theorem. 

In §3, we reduce the Euler system into a second-order nonlinear elliptic equation 
in the subsonic region and then reformulate the wedge problem into a one-phase 
free boundary problem for the second-order nonlinear elliptic equation with the 
shock-front as a free boundary. 

In §4, we use the hodograph transformation to make the free boundary into a 
fixed boundary, in order to reduce the difficulty of the free boundary. After that, 
we only need to solve for the unknown entropy function A as a fixed point. 

In §5, for a given entropy function A, we solve the reformed fixed boundary 
value problem in the unbounded domain and determine a priori the downstream 
asymptotic function of the potential function at infinity. Then, in §6, we prove that 
the entropy function A is a fixed point via the implicit function theorem, which is 
one of the novel ingredients in this paper. 

In §7, we determine the decay of the solution to the asymptotic state in the 
physical coordinates. 
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In §8, we establish the stability of the transonic solutions and transonic shocks 
under small perturbations of the incoming flows and wedge boundaries. We finally 
give some remarks for the problem when the downstream state of the background 
solution is on arc TH in §9. In Appendices, we show two comparison principles 
and derive a criterion for different arcs TS and TH on the shock polar, which are 
employed in the earlier sections. 

Finally, we remark that the stability of conical shock-fronts in three-dimensional 
flow has also been studied in the recent years. The stability of conical supersonic- 
supersonic shock-fronts has been studied in Liu-Lien [24] in the class of BV solutions 
when the cone vertex angle is small, and Chen [8] and Chen-Xin-Yin [12] in the class 
of smooth solutions away from the conical shock-front when the perturbed cone 
is sufficiently close to the straight-sided cone. The stability of three-dimensional 
conical transonic shock-fronts in potential flow has been established in Chen-Fang 
[2] with respect to the conical perturbation of the cone boundary and the upstream 
flow in appropriate function spaces. Also see Chen-Feldman [3]. 

2. Mathematical Setup and the Main Theorem. In this section, we first 
formulate the wedge problem into a free boundary problem for the composite-mixed 
Euler equations, and state the main theorem. 

As is well-known, for a uniform horizontal incoming flow Uq = {u^q,0,Pq , p^) 
past a straight wedge with half-wedge angle Oq, the downstream constant flow can 
be determined by the Rankine-Hugoniot conditions, that is, the shock polar (see 
Appendix B and Fig. 2.1). According to the shock polar, the two flow angles 
are important: One is the critical angle 9^ that ensures the existence of the at¬ 
tached shocks at the wedge vertex, and the other is the sonic angle 9% for which the 
downstream fluid velocity at the sonic speed in the direction. When the straight 
wedge angle 9^ is between 0^ and 0^, there are two subsonic solutions; while the 
wedge angle 9^ is smaller than 0^, there are one subsonic solution and one super¬ 
sonic solution. We focus on the subsonic constant state Uq = (uq ,pg ,Pq) with 
Uq • (sin00, — cos00) = 0. Then the transonic shock-front Sq is also straight, de¬ 
scribed by xi = soa: 2 . The question is whether the transonic shock solution is stable 
under a perturbation of the incoming supersonic flow and the wedge boundary. 

Assume that the perturbed incoming flow U~ is close to Uq , which is supersonic 
and almost horizontal, and the wedge is close to a straight wedge. Then, for any 
suitable wedge angle (smaller than a critical angle), it is expected that there should 
be a shock-front which is attached to the wedge vertex. If the subsonicity condition 
is imposed in the far held downstream after the shock-front, then the flow U between 
the shock-front and the wedge should be subsonic. Since the upper and lower 
subsonic regions do not interact with each other, it suffices to study the upper part. 

We now use a function b{xi) to describe the wedge boundary: 


dW = {x G : a ;2 = b{xi), 0(0) = 0}. 


( 2 . 1 ) 


Along the wedge boundary cIW, the slip condition is naturally prescribed: 



( 2 . 2 ) 


Let the shock-front 5 be Si = 17 ( 0 : 2 ) with 17 ( 0 ) = 0. Then the domain for the 
subsonic how is denoted by 


fls = {x G : xi > ( 7 ( 0 : 2 ), X 2 > 0 ( 0 : 1 )}, 


(2.3) 
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sonic .circle 


o 



Figure 2.1. Two arcs TS and TH on the shock polar 


and the shock-front S becomes a free boundary connecting the subsonic flow (ellip¬ 
tic) with the supersonic flow (hyperbolic). 

To be a weak solution of the Euler equations (1.1), the Rankine-Hugoniot con¬ 
ditions should be satisfied along the shock-front: 


[pUl ] = (7'{X2)[PU2], 

[pul +p]= a'{x2)[puiU2], 

[pUlU2] = Ct'{x2)[pU2'^ + p ], 
[pu,{E+J^^)]=a'{x2)[pu2{E+I^^)] 


(2.4) 




as the free boundary conditions on S, where [ • ] denotes the jump between the 
quantity of two states across the shock-front. 

For a fixed uniform supersonic state Uq , there is an arc on the shock polar 
corresponding to the subsonic states. When the wedge angle is less than the critical 
angle 9^ > the tangential point T corresponding to the critical angle divides 
arc EIS into the two open arcs TS and TEI. The nature of these two cases is very 
different. 

In this paper, we analyze the existence, stability, and asymptotic behavior of 
steady transonic flows with a transonic shock in the important regime TS for the 
wedge angle 9^. To state our results, we need the following weighed Holder norms: 
For any x, x' in a two-dimensional domain E and for a subset P of dE, define 


(5x := min {dist(x, P), 1} , 


mm{(5x,^x'}, 

= mm{|x| -I- 1, |x'| -I- 1} 
= minjAx, Ax'}- 


Ax := |x| -h 1, 

Ax := dist(x, P) -|- 1, 


Ax.x' 


Ax.x' 
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Let a G (0,1), cr,T, I G M, and let fc be a nonnegative integer. Let k = (^ 1 ,^ 2 ) be an 
integer-valued vector, where fci, ^2 > 0, |k| = ki + ^ 2 , and We define 


[fto-iriYE ■■= sup (<5r^{'=+'^'°>A;AL+'=|i^V(x)|), 

|k|-fc 

[/]fc.a;(r.0;E - ^Up 
x.x GE 
x^x' ,\k.l=k 


cmax{fc+a+(T,0} * r Al+k+a l-^^f (^) 
x,x' x,x' x,x' |x-x'|“ 


k 

11*^ ''k,of,{T,l)',E ' / j I-'' \i,0‘,{T,l)‘,E ' !-*' -lfc,a;(r,Z);£" 

2=0 


(2.5) 


For a vector-valued function f = {fi, f 2 ,''' 1 fn), we define 


\{W,P 

\k,a-(T,l)\E 


n 


:=Ell/.ll 


L^):P 

k,a;(T,l);E' 


2=1 


For a function of one variable defined on (0, 00 ), we define the Holder norms with 
a weight at inhnity. The definition above can be reduced to one-dimensional if we 
keep only the weights at infinity. Then the notation becomes ||/||fc,c(;(T);(o,oo)- 
We also need the norms with weights at infinity which apply only for the deriva¬ 


tives: 


lk,a;(T,l);E 


■■= \\f\\cOiE) + \\Df\\ 


{<y+l)-,P 

fc—l,Q:;(r,Z + l);£^' 


( 2 . 6 ) 


Similarly, the Holder norms for a function of one variable on (0,oo) with only the 
weights at infinity are denoted by ||/||fe,„;(r);(o.oo)- 

In terms of supersonic flows, we prescribe the initial data: 


C^|i = Uo{x 2 ) on I := {a;i = 0}. 
Let 12“ be the domain for the incoming flows defined by 


12“ = {x : 0 < Si < 250 X 2 } . 

(2.7) 

For a given shock S = {xi = cr(x 2 )}, let 


12^ = {x : 0 < xi < (j{x 2 )} ■ 

(2.8) 


We now fix parameters a, P G (0,1) with suitably small 13, depending on the back¬ 
ground states. 

Then we can conclude that there is e > 0, depending on the background states, 
such that, when 


l|C ^0 - C^o Il2,a;(i-L/3);I < £ for some ^ > 0, (2.9) 

there exists a constant Cq > 0, independent of e, and a unique supersonic solution 
U~ = [vT ,p~, p~){x,y) of system (1.1) with the initial condition U~\x = Cfo, well 
defined on 12“, such that 


\\U- - C/o-||2.c;(l+/3);a- < - C/o" ||2..;(1+/3);I. (2.10) 

This can be achieved by rewriting the problem as an initial-boundary value prob¬ 
lem in the polar coordinates [r,6) so that system (1.1) is still a hyperbolic system, 
domain 12“ becomes a half strip with 6 time-like and r space-like, the initial data 
is on {r > 0, 0 = 0}, and the boundary data u = 0 is on the characteristic line 
{r = 0,0 < 0 < arctan(2so)}. This is a standard initial-boundary value problem 
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whose almost-global existence of solutions can be obtained as long as e is sufficiently 
small. 

Assume that the wedge boundary satisfies 


11^ ~ ^o||nQ;(/3);R+ < £• 


( 2 . 11 ) 


Theorem 2.1 (Main Theorem). Let the background solution {Uq , Uq } satisfy that 
Uq is on arc TS in Fig. 2.1 for the straight wedge boundary X 2 = bo(xi) = 
tan0o3ii,xi > 0. Then there is e > 0 such that, when the initial data Uq and 
the wedge boundary dW = {x 2 = b{xi),b{0) = 0} satisfy (2.9) and (2.11) re¬ 
spectively, there exist a strong transonic shock S := {a;i = o'(x 2 )}, a transonic 
solution {U~,U} of the Euler equations (1.1) in Lis, OLnd an asymptotic down¬ 
stream state U°° = (u°°,p(j", p°°) = V°°{x 2 — tanSg 2 ^ 1 ) for an appropriate function 
V°° : [—e, 00 ) —>■ with • (sin^o, — cos 6 * 0 ) = 0 for the wedge angle 9o such that 

(i) U~ is a supersonic flow in Ll^, and U is a subsonic solution in Lis; 

(ii) The Rankine-Hugoniot conditions (2.4) hold along the shock-front S; 

(iii) The slip condition (2.2) holds along the wedge boundary dVL; 

(iv) The following estimates hold: 


\\U- - + W " II + 11^ - Un\o,.;ip,iy,ns 

+ ||(t'(-) - - C^o'llo,a;(l+/3);[-e,oo) 


< C (||17o - Uq ||2,a;(l+/3,0);I + 11^ “ ^0 II l,a;(/3);R+) > 


( 2 . 12 ) 


where C is a eonstant depending only on U^, but independent of e. 

Moreover, solution U is unique within the class of transonic solutions such that 
the left-hand side of estimate (2.12) is less than Ce. 


Remark 1. Estimate (2.12) implies that the downstream flow and the transonic 
shock-front are close to the background transonic solution. The subsonic solution 
U converges to U°° at rate |x|“^ and the slope of the shock converges to the slope 
of the background shock at rate 


Remark 2. Theorem 2.1 indicates that the asymptotic downstream state U°° gen¬ 
erally is not a uniform constant state. If the x-coordinates are rotated with an¬ 
gle 6q into the new coordinates {xi,X 2 ) so that the unperturbed wedge boundary 
i9>Vo = {x 2 — tan 6 * 0 X 1 } becomes the xi-axis: 

(xi, X 2 ) = (cos 6*0 xi -I- sin 6*0 X 2 , — sin 6*0 xi -f cos^ 0 X 2 ), 

then = V°°{x 2 ). In Lagrangian coordinates, y = {yi,y 2 ), determined by 
(3.1) in §3, the asymptotic downstream state is a function of y 2 in general: U°° = 
U°°{y 2 ) = (u°°( 2 / 2 ),Po’, P°°( 2 / 2 ))- However, our argument also shows that, in the 
isentropic case with a constant Bernoulli quantity B (see (3.14)), the asymptotic 
state must be uniform and equal to the background state. Also see Chen-Chen- 
Feldman [1]. 


3. Reduction of the Euler System and Reformulation of the Wedge Prob¬ 
lem. In this section, we hrst reduce the Euler system into a second-order nonlinear 
elliptic equation and then reformulate the wedge problem into a free boundary prob¬ 
lem for the nonlinear elliptic equation with the shock-front as the free boundary. 
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From the first equation in (1.1), there exists a unique stream function -0 in domain 
U such that 

= {-pu2,pui) 

with '0(0) = 0. 

To simplify the analysis, we employ the following coordinate transformation to 
the Lagrangian coordinates: 


f 2/1 = xi, 

1 2/2 = '0(a;i,X2), 


(3.1) 


under which the original curved streamlines become straight. In the new coordinates 
y = (2/1,2/2), we still denote the unknown variables [/(x(y)) by U{y) for simplicity 
of notation. 

The original Euler equations in (1.1) become the following equations in divergence 
form: 


(—) -(-) =0^ 

^pui^y^ 

(3.2) 

pui^y^ ^ ui ^y^ 

(3.3) 

{U2)y, +Py 2 = 0, 

(3.4) 


(3.5) 


Let T : 2/1 = < 5 '( 2 / 2 ) be a shock-front. Then, from the above equations, we can 
derive the Rankine-Hugoniot conditions along T : 


[—] = -p]a'(j„), 

(3.6) 

h + —] = -[—] ^'(2/2), 

'■ pui ■' Ui 

(3.7) 

[U2] = [p]ct( 2 / 2 ), 

(3.8) 


(3.9) 

The background shock-front now is To : 2/1 = S12/2, with = 

= Poutoij^-t'Aneo) > 

0. Without loss of generality, we assume that the supersonic 
domain D“ defined by 

solution U~ exists in 

D“ = {y : 0 < 2/1 < 2 si 2 / 2 } . 

(3.10) 

For a given shock function ^{1/2), let 


= {y : 0 < 2/1 < o'( 2 / 2 )} , 

(3.11) 

Da = {y ; 0 < 2/2, <7(2/2) < 2/1} ■ 

(3.12) 


In either the supersonic or subsonic region, X 2 can be solved as a function of y 
since ipx 2 = P^i ^ 0. Let X 2 ■= 0(y) in the subsonic region and X 2 ■= 0~(y) 
in the supersonic region Dr. Given U~, we can find the corresponding function 
4>~. We now use function 0 (y) to reduce the original Euler system to an elliptic 
equation in the subsonic region. 
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By the definition of coordinate transformation (3.1), we have 

Yyi — 5 ^ 2/2 ~ ? 

Ui pUl 

that is, (j){y) is the potential function of the vector field (^, ^)- 
Equation (3.5) implies Bernoulli’s law: 



IP 

( 7 - 1)P 


5 ( 2 / 2 ), 


(3.13) 


(3.14) 


where B = B{y 2 ) is completely determined by the incoming flow U~ at the initial 
position I, because of the Rankine-Hugoniot condition (3.9). 

From equations (3.2)-(3.5), we find 


(blnp- lnp)yi = 0, 


which implies 


p = A{y 2 )p'^ in the subsonic region Do-- 


(3.15) 


With equations (3.13) and (3.15), we can rewrite Bernoulli’s law into the following 
form: 

_|_ 1 ry 


+ = Bp^. 


20 L ' 7-1 


(3.16) 


In the subsonic region, |u| < c := Therefore, Bernoulli’s law (3.14) implies 


p^ 1 > 


2 ( 7 -1)5 


(3.17) 


7 ( 7 + 1 )^' 

Condition (3.17) guarantees that p can be solved from (3.16) as a smooth function 
of {A,B,V(j)). 


Assume that A = A{y 2 ) has been known. Then {u,p,p) can be expressed as 
functions of Vcj): 


P = p{A,B,\/(j)), 


u=( 


1 4 >yi 

P^’yi P4’y2 


), 


p = Ap'f, 


(3.18) 


since B = 11 ( 2 / 2 ) is given by the incoming flow. 

Similarly, in the supersonic region D“, we employ the corresponding variables 
(A“, B, (j)~) to replace 17“, where B is the same as in the subsonic region because 
of the Rankine-Hugoniot condition (3.9). 

We now choose (3.4) to derive a second-order nonlinear elliptic equation for cj) so 
that the full Euler system is reduced to this equation in the subsonic region. Set 

N^=u 2, N^=p. (3.19) 

Then we obtain the second-order nonlinear equation for cf: 

{N\,+{N%,=0, (3.20) 

where 77* = 77*(A(2/2), 5 ( 2 / 2 ), V^), i = 1,2, are given by 


N\A,B,V4>) 

N\A,B,V4>) 


__ 

</’y2P(^(y2), 5(2/2), v<)))’ 

^( 2/2 )p(A( 2/2 ), 5 ( 2 / 2 ), V(())^. 


(3.21) 
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Let q = \/ul+ U 2 - Then a careful calculation shows that 


1 _ui{c^-uj) 

’ 

(3.22) 

Nl =Nl = 

^y2 ^vi 

(3.23) 

2 2 2 

2 _c p q ui 

(.2 _ g2 ■ 

(3.24) 

_2 _2. 2 

Nl Nl = / i >0 

‘Py2 ^y2 ‘^yi 

(3.25) 


Thus, the discriminant 


in the subsonic region with pui ^ 0. Therefore, when solution (jj is sufficiently 
close to (/>o (determined by the subsonic background state Uq) in the norm, 
equation (3.20) is uniformly elliptic, and the Euler system (3.2)-(3.5) is reduced to 
the elliptic equation (3.20) in domain D^., where a is the function for the transonic 
shock. 

The boundary condition for (j) on the wedge boundary {y 2 = 0} is 


<(>( 2 / 1 , 0 ) = h{yi). 


(3.26) 


The condition on T is derived from the Rankine-Hugoniot conditions (3.6)-(3.8). 
Condition (3.6) is equivalent to the continuity of (j) across T": 


It also gives 


[(j)]\j- = 0. 


<3-'(2/2) = -|^(ct(2/2),2/2)- 


(3.27) 

(3.28) 


Replacing (j'{y 2 ) in (3.7) and (3.8) with (3.28) gives rise to the conditions on T: 

1 


G(t/-,AV0) = 


P4‘y2 

H{U-,A,V^) = + [K][N^] = 0. 


- [<(>y2][^P'^?^'yi] = 0, (3.29) 

(3.30) 


We will combine the above two conditions into the boundary condition for (3.20) 
by eliminating A. 


By calculation, we have 


/Vl - ^ P'* ^^2 

■^'■‘4 1 2 2 ’ 

-i - Ic^ - 


Nl = - 




2 2 


(3.31) 

(3.32) 


Thus, we obtain 


Ha — H\[(l)y^] + N'\[(l)y^] 


1 

eT 

to 

'U2' 

-Ul- 


1 

7 - 1 c2 - 

(? — q^ 

.pu\. 


> 0 , 
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and 


Ga 


< 


[K](^ + ^V^Na) - [4’y2]KNA 

\ (Py^ / 


U2p''{q^ + ^P 


Ui{c^ — q‘^) lpui\ Ui{c^ — q^) 


o7-l 


(u2'^ + 



.Ml. 


since 1 < 0 and M 2 _ is close to 0. 

Therefore, both equations (3.29) and (3.30) can be solved for A to obtain A = 
gi{U ~and A = g 2 {U~,\7(j)), respectively. Then we obtain our desired condi¬ 
tion on the free boundary {Le. the shock-front): 


g{U-y<l)) := (32 - 5 i)(C^-,V</>) =0. (3.33) 


Then the original transonic problem is reduced to the elliptic equation (3.20) with 
the fixed boundary condition (3.26) and the free boundary conditions (3.27) and 
(3.33), and A is determined through either of (3.29)-(3.30). 


4. Hodograph Transformation and Fixed Boundary Value Problem. In 

order to reduce the difficulty of the free boundary, we employ the hodograph trans¬ 
formation to make the shock-front into a fixed boundary. After that, we only need 
to solve for the unknown function A. 

We now extend the domain of <()“ from D“ to the first quadrant D“ U Do-. Let 
(/>Q = _ 32 ; which is the background potential function. We can extend 4>~ into 

Pq '^20 

D” U Da- such that 


(f) = when 0 < 2 si ?/2 < 3i - 1- 

We then use the following partial hodograph transformation: 

Zl =(j)-((>-, 

Z2 = 2/2, 

SO that 2/1 is a function of {zi,Z2)' yi = (p{zi,Z2)- 
Let 

M\U-,A,W(I)) = N\A,B,V(I)) +N'^{A,B,V(j>) 


[^^ 2 ] 

[4’yi] 


and 


M^{U-,A,V(j)) = 


M\z,A,(p,W(p) 


N^{A,B,\J(j)) 


[<A: 


yi\ 


(4.1) 


=-M\U {g},Z2),A,dy^4> {ip,Z2) + —,dy^(l) {ip,Z2) - i = l,2. 

Therefore, equation (3.20) becomes 

(M^{z,A,ip,Vip)) + (M^(z,A,ip,^ip)) = 0 . 


(4.2) 













Notice that 
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Also 




K., = 


^Nl^+2NlJc^y,][4>y,] 

+ [4’yp ^ 


’ry2 

Vt^yi] 

+ Nljllly2\+N\ 


K„ = 

Nl 

'ry2 


+ Nlj<t>y,]-N\ 



Nl 

'ry2 




1 /_ _1 _ _2 \ 


„-(N„fl>o. 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


which implies that equation (4.2) is uniformly elliptic, for any solution ip that is 
close to ip^ (determined by the background solution Uq ) in the norm. 

Then the unknown shock-front 'T becomes a fixed boundary, which is the 2 : 2 -axis. 
Along the Z 2 -axis, condition (3.33) is now 

g{z,ip,Vip) = g{U~{(p, Z 2 ),dy^(l)~{(p, Z 2 ) + — ,dy^(j)~ {(p, Z 2 ) - —) = 0. (4.7) 

<Pzi ‘fzi 

We also convert condition (3.30) into the z-coordinates: 

H{z,A,ip, Vip) 

:= H{U~ {(p, Z 2 ), A, dy^cl)~ {ip, Z 2 ) + —,dy^(l)~{ip,Z2) - —) = 0 (4.8) 
along the 2 : 2 -axis. 


The condition on the Zi-axis can be derived from (3.26) as follows: Restricted 
on Z 2 = 0, the coordinate transformation (4.1) becomes 


= Kyi) - ^-(l/uO)- 


Then yi can be solved in terms of zi so that 

yi = ^{zi,A) = h{zi). (4.9) 


Let Q be the first quadrant. Then the original wedge problem is now reduced to 
both solving equation (4.2) for ip in the unbounded domain Q with the boundary 
conditions (4.7) and (4.9) and determining A via (4.8). 


This will be achieved by the following fixed point arguments. Consider a Banach 
space: 

X = {A : A(0) = 0, < 00 } 

as defined in (6.3) below. Then we define our iteration map JT '■ X — X through 
the following two steps: 

1 . Consider any A = A{z 2 ) so that A — wt £ X satisfying 

( 4 . 10 ) 

for some fixed constant Cq > 0, where Wt = Wt{z 2 ) is determined by (6.2) below. 
With this A, we solve equation (4.2) for ip = ipA in the unbounded domain Q with 
the boundary conditions (4.7) and (4.9) in a compact and convex set: 

Ss = {‘P ■ Wp — P°°\\ 2 ^aK,p^o)^ — sufficiently small i5 > 0 (4.11) 
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in the Banach space: 

B={^p : \\tf- 0 < a' < a, 0 < /3' < ^, (4.12) 

where ip is determined by (5.2). Equation (4.2) is uniformly elliptic for p G Tis for 
small ^ > 0. The existence of solution pA G Es will be established by the Schauder 
hxed point theorem in §5. 

2 . With this p = pA, we solve (4.8) to obtain a unique A that defines JT(A—Wt) = 
A-Wf 

Finally, by the implicit function theorem, we prove that has a fixed point 
A — wt in § 6 , for which A satisfies (4.10). 


5. An Elliptic Problem to Determine p in Domain Q. In this section, for 
given A satisfying (4.10), we solve equation (4.2) for p in the unbounded domain 
Q with boundary conditions (4.7) and (4.9). Before this, we determine a priori the 
limit function p^ at infinity. 


5.1. Determine a priori the limit function p°° at infinity. First, we assume 
that the asymptotic downstream state depends only on y 2 , which will be verified 
later. Then we determine the limit function (jf° for (p. From (3.2), we expect the 
flow direction at infinity is the same as that of the wedge. That is, 

"“2 ^ . uto 

(py^ = - > tant^o = -rr yi ^ oo. 

U\ ^10 


Then 


(p°° = tan^oyi + l{y2)- 

Replacing (p with p°° in Bernoulli’s law (3.16), we obtain 


(tan^o)^ + 1 


-J—Ap'-f+^ = Bp^. 

7-1 


From (3.4), we expect that pressure p ^ Pq and then relation (3.15) becomes 
Pq = A{p°°)^ so that A = A{y 2 ) and p°°{y 2 ) = ( Afas ) ■ Therefore, the above 
equation becomes 


(tan^o)^ + 1 

2iny2)f 


7 

7 - 



(7+l)/7 



2/7 


(5.1) 


This equation gives the expression for l'{y 2 )- We can find /(j/ 2 ) by integration with 
Z( 0 ) = Wo, where wq is the limit of 6 — 60 as —>■ 00 . 

Then we employ 

= -p-){p°°,Z2) (5.2) 

to solve for p°°. Also, equation (5.2) restricted on Z 2 = 0 gives rise to 


Zi=tan0obo + wo-p ( 60 , 0 ), 


from which we can solve for bo- 

By the definition of p°°, we know that p°° satisfies (4.2). That is. 


(5.3) 
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5.2. Linearization. Let 

^5 = {w : 

where the wedge boundary dW is the zi-axis. We will omit dW in the norm when 
no confusion arises. 

To solve equation (4.2) in the first quadrant Q, we first linearize (4.2) and solve 
the linearized equation in bounded domains, and then take the limit to obtain a 
solution in the unbounded domain Q. 

For given ip such that p — we define a map Q: p — p°° =: Q{p — p°°) 

in Tjs as the solution of the linearized equation (5.5) of the nonlinear equation (4.2) 
at p as described below, and show that there exists a fixed point that is a solution 
for equation (4.2). 

We use a straight line := {z 2 = R— zi} to cut off Q into a triangular domain 

:= {0 < Z2 < R - zi,zi > 0}. 

Let 

v = p-p°°, ( = b-bo. 

Taking the difference of equations (4.2) and (5.3) and linearizing the resulting equa¬ 
tion lead to 

Lt=1.2 

where 

af^ = . (z, A, + s(<^ - ^°°), V(^“ + s(^ - p°^))) ds, (5.6) 

M;(z, a + s((^ - ^-), + s(<p - P^)) ds (5.7) 

Jo 

for ij = 1,2. 

Replace p by p°° in af^ and denote by a“. Then afj — afj is bounded in the 
Holder norm || • |li is bounded in || • 111 ^Iso, the uniform 

ellipticity of equation (5.5) follows from (5.6) and the uniform ellipticity of (4.2) for 
the solutions close to p^, provided that S in (5.4) is chosen sufficiently small. 

The boundary condition on the zi-axis is 

vU=o = C (5.8) 

On the cutoff line we prescribe the condition; 

vIlr = C(i?), (5.9) 

which is compatible with the condition on the zi-axis at point (i?,0). 

Condition (4.7) on the 2 : 2 -axis can be linearized as follows: Condition (4.7) can 
be rewritten as 


g{z,p,\/p) -g{z,p°°,yp°°) = -g{z,p°°,yp°°). 
Therefore, we derive the oblique condition: 

X] =: go, 

i^l,2 


(5.10) 
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where 


-^o 

+s(<P - V((/7°° + s(<p-<p“)))ds, 

Jo 

which have all the corresponding bounded Holder norms. 

When Uq is on arc TS, the direction of !>< = {vi,i' 2 ) is 


VI = 9v., 




(7 - 1)uIGaHa{c^ - q^) ^^pui 


( - 2c^pUlM2[——] + U2(c^ - ^ 1 )) 

V'-pui'' pwi J 


> 0 

since [-^\ < 0, and 


V2 = = 


-p'^ ^U2 


(7 - 1)uiGaHa(.c^ - q^) 


_ n 

2 _ „2'l 


where 


Cp = [p] (c^ + (7 - l)q^ - 7 U?) + (7 - 1)P9^M2 + [ - ] P^c^Uiq^. (5.11) 

^ ' puG 

Since, on arc TS, Gp < 0 from (B.ll), we have 


V2 < 0 . 


In particular, if S is small, then 

1 , 

1^2 < 2*^20 < 0^ 

where i'^q is the quantity 122 for the background subsonic state. This implies that 
condition (5.10) is uniformly oblique. 

Set 

e = l|t^0 - Uq ||2,0!;(1+/3);D- + \\b - &0 || 2,O!;(,S);(0,oo)' 

Now, for any function / of {U~,U), we use / to denote the value at the back¬ 
ground states: / = f{UG ,Uq). We also omit domain and boundary dW in the 
norms when no confusion arises. In the following estimates, all generic constants 
C’s and cq are only dependent on the background states and a, /3, but independent 
of e, S. 


5.3. C° estimate for v. We employ the comparison principles. Theorem A.l and 
Theorem A.2, to estimate v. 

We decompose matrix A = (Sy) into A = KK^, where 

( / aiia22-g^„ ai2 \ 

V 022 1 , 

0 \/a22 J 

and = 1,2, are the coefficients of the second order terms in the linearized 

equation of the nonlinear equation (4.2) at the background solution. 
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We define the transformation z = Kz, where z = ( 21 , 2 : 2 ) is a new coordinate 
system. Then X)i=i 2 ^ 0 % = implies 

Li=1.2 

We use the polar coordinates (r, 6) for z to construct a comparison function for 
V. That is, 

r = |z|, 0 = arctan(—). 

Let 9 = t9 + T. Define 

V = sin0, (5-12) 

where r > 0, and s and t will be chosen later. 

We compute 

= A^v = (s^ — t^) sin0. (5.13) 

ij = 1.2 

Let s = —/3, t = a, and 0 < ,5 < a in (5.12). We set vi = sin(a0 + r). 

First, (4.10) leads to 

A'{Z 2 ) = 0{e) max( 22 ,1)“^“^ min( 22 ,1)““^ 

= 0{e) max( 22 , l)“““^“^[max( 22 ,1) min(z 2 ,1)]““^ 

= 0(e) max(22,1)“““^“^Z2 
= 0{e){z2 + l)-“-^-V“-ysin6l)“-\ 


Since (4.11) gives \\ip — II 2 a-(/ 3 ^oy^ < ^ , by the expressions of afjyb'j (see 
(5.6) (5.7) ), we have 

\afj - tty I < CS. 

(<■).. = O(<5)r-'5-2(sin0)“-i + 0(1)A'(22) 

= 0(<5)r-^-2(sin6»)“-i + 0(£)(z2 + 1)-“-^" V“-ysin6l)“-b 
bf = 0{e)r-^-^. 

= 0{e)r-P-^ + 0{e)r-P-^A'{z2) 

= 0{e)r~^~^ + O(e)r“^“^(sin0)““b 


(a 




- = 0{S)r 

{afj)z^d^,vi = 0{5)r 
< 0{5)r 
= 0{5)r 
b^dziVi = 0{e)r 
{bf)z,vi = 0{e)r 


- sin 9. 

-P-2^ 

-2/3-3(sin6l)“-i + 0{e){z2 + 

-'^-2(sin6l)“-i + 0(£:)22"“r“-^-ysin6»)“-i 

“^“^(sin6l)““^ + O(e)r“^“^(sin0)“^. 
-2/3-3 

-'5-ysin6»)“-b 


(sin6»)“-i 
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Therefore, we compute 

L'^vi = + iafj - aij)dl,.)vi + + b^d^.Vi + {bf)^,Vi 

< —cor~^~^ + 0(s) r~^~^(smd)~^, 


where cq > 0 is a constant depending on the background state and a, /3. 

Let V 2 = r~P sin“ 9. In the same manner, we can obtain the following estimate: 


L‘^V2 < —cir 


-'5-2sin“-2 0. 


Set Vs = vi +V 2 - Thus, we have 

L'^vs < —cor~^~^ < 0. 


By Theorem A.l, achieves its positive maximum on the boundary. 
On Z 2 = 0 and — < Ce. 

Let 6*0 = arctan(—1^). We compute VzU on the 2;2-axis: 

Ujj = ^ (s cos 00 sin 0 — t sin 0o cos 0), 

= r^~^{s sin 0o sin 0 + t cos 0o cos 0). 


Then 


VzU = (K -^)' VzU = r 


s COS ^0 sin 9—t sin 9q cos 0 
„s-l I fell 

s sin 9 
k 22 sin 9o 


, , /3cos0o sin0 + Q!sin0o COS0 + (a + /3) sin“ 00 COS0O 

(wajzi =-^- r 


- 0-1 


I'll 


^ /3(sin0 + sm" 0o) .-i 

^ fc22Sin0o 

Then, when /3 is suitably small, we have 

Di^ivi + V2) < 


Assume that ^ achieves its maximum eM at some point P on the Z2-axis. We 
know that Di^{^){P) < 0. 

Since |go| = \g{z,(po,'^‘Po) — ff(z,< Cer~^~^, we obtain that, at 
point P, 

0 > DuV - —D^vs 

V3 

= go — c'^V — £MDi,Vs 

> —C£{l + £M)r~^~^+M£Cir~^~^. (5-14) 

This implies that M < ^ for sufficiently small e. 

Thus, we obtain the estimate for v. 

|u| < C£r~^. 


(5.15) 
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5.4. C^'°‘ estimate for v at corner O. In (5.12), let s = 1 + a and t = 1 + a + r. 
We define 

114 = sin ((1 + a + t)0 + r). (5.16) 

By (5.13), it is easy to check 

< -C2r“-i. 


On the Z2-axis, we have 

D^va = r“ ("^((a + 1) sin((Q; + t)0o + r) - r sin0o COS0) + 

V fell ^ fc22 sin 6*0 / 

< -cr“, 


provided that a and r are suitably small. 

Then we can use eC'i ;4 as a comparison function to control w = v—v{Q)—Dzv{0)-'z 
for r < 2. 

Denote any quarter ball Br(0) 0 Q with radius r by In B^, 

L‘f’w = '■=- (X! + D^viO) ■ z))^. 

j^l,2 i^l,2 2^1,2 

> -C'er“-i > L'^iCevi). (5.17) 


By Theorem A.2, we have 


sup 




+ 1 . 


On dB 2 n {{z 2 = 0} U {|z| = 2}), we see that < C. 

Assume that achieves its maximum M at a point P on the 22 -axis. The 
oblique condition (5.10) implies 

"Y = go = 0(er“). 

i=l,2 


The same argument 


as in (5.14) implies that, at the maximum point P, 


0 > D.(-) 

Vi 


1 

Vi 


[Du'w - DuVa] 

^ ?;4 ^ 


> —( —c'^w — eCr°‘+ eMcor°‘), 

Va 

which implies that M < ^. Thus, w < eCr^'^°‘ in 
Similarly, we obtain the corresponding lower bound. 
Therefore, we conclude 


|?i'(z)| < eCr^~^°‘ for any z € B^ ■ (5.18) 

With estimate (5.18), we can use the scaling technique to obtain the C^’°‘ estimate 
for w up to the corner. More precisely, for any point P* € B^ with polar coordinates 
(d«, 0,), we consider two cases for different values of 0,. 


Case 1: Interior estimate for 9^ € Set Pi = B^{P^) and P 2 =B^{PC)- 

6 3 

Then Pi C P 2 C P^. By the Schauder interior estimates (c/. (4.45) and Theorem 
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8.33 in [17]), we have 

1 = 1,2 

where fi is defined in (5.17), C is a constant independent of d*, and the weight of 
the norm is up to dB 2 . Therefore, by (5.18), we conclude 

||w||i,„;Si < < Ce. (5.19) 

Case 2: BouTiddTy cstiuidtG fov 6^ ^ ^ ot 0^ <i ^. Denote = Q D B 2 d^ (p*)- 
By the Schauder boundary estimate (c/. (4.46) and Theorem 8.33 in [17]), we have 

IL„l|(0) 


I 


< c(||»||„,„,* + y; ||/,||‘«,B. + + lISollJl, 

2 ^1,2 

< eCd]:+“. 

Combining Case 1 with Case 2 yields the corner estimate: 

ll^lll.a;S+ = Ik + ^(0) + ■ z|li,„.B+ < Ce. 


B3n{zi=o} 


(5.20) 


5.5. C^’°‘ estimate for v away from corner O. Away from the corner, suppose 
z* = {zl,z^) € with polar coordinates {R*,0*) for 1 < i?* < f. Two cases will 
be considered below. 

Case 1. Z 2 > 5 - Let i?, := (z*) fl Q. If 6* > arctan2, we employ the 

~T 

Schauder boundary estimate; for 9* < arctan2, we employ the Schauder interior 
estimate. Notice that 

K/oI.B, < (4)“[<]o.c;H. < Ce < 1. 

Therefore, where A > 0 is a constant only depending on the back¬ 

ground states. Hence the constant C in the following estimate is independent of 


kllS;H. < C(1|U1|0,0;H. + ||50||^'\.. 3 .. 


< C'e(i7*)-^, 

where the weight in the superscript is the distance up to {dB^* (z*))nQ. Shrinking 

~T 

the domain by setting H' := B (z*) n Q yields 

T 

kl|l,a;(,a, 0 );B' < Cs. 

Case 2. 2 :| < 5 . In this region, we apply Schauder boundary estimate in B^ := 
Bi{z*) Pi Q to obtain 

kllLi;B+ ^ c(||u||o,o,3+ + ||C||l.„;(z-l.z* + l)) 

< Ce(i?*)-^, 


I^lll.a;(/3,0);B+ — 


which implies that 
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Combining the estimates in the two cases above, we obtain the estimate for v in 

lklll,a;(/3,0);Q^/2 ^ (5-21) 

5.6. C^’“ regularity. For the (7^’“ estimates with a weight to the 2 ;i-axis, we 
rewrite equation (5.5) into the non-divergence form: 

+ {bt)ziv) = 0, (5.22) 

ui=i.2 

with the boundary condition on the 2 : 2 -axis: 

Y (5-23) 

where 

gi = go- c‘^v + Y^^^~ ■ 

i=l,2 

We use the same scaling as in the (7^’“ estimates and consider two cases. One 
is away from the 2 i-axis and the other is close to the zi-axis. We will estimate the 
(72’“ norm in Q^/^. Suppose z* = {z*,z*) G 

Case 1. Z 2 > 1. S* := i?** (z*) O (J is the same as in the (7^’“ estimates. Set 

~T 

T = f|{zi = 0}. 

Similar to the situation in Section 5.5, because the coefficients in (5.22) have the 
proper decay rate in Z 2 direction, we have 

The Schauder interior and boundary estimates (c/. Theorem 6.26 in [17]) imply 

<C'(||u||o.O;R. + ||5l|ll!i;T)- (5-24) 

Then 

||5l||£;T < C'(l5o|^^;T+e|MI?iB.) 

leads to 

lkll2,O!;(,fl,0);Bi < Cs. (5.25) 

Case 2. z* < 1. Define '!;(z) = v(z) — v(zi,0) — Vv(z*,0) ■ z and let = 
Bk(z*) n Q. Then 

||^)||0,0;B. < 

< Ce\z*\-f^zY°‘- 

The same Schauder interior and boundary estimates as (5.24) applied to v gives rise 
to 

II«I|2.c.;(o.o);b: <c'c|z*|-^2;^+“. 

Hence, we conclude that 

^ C-c. (5.26) 

Estimate (5.25),(5.26), together with (5.15), leads to 


IMI 


(-l-a);aw 


2.a;(/3,0);Q«/^ — 


< Ce. 


(5.27) 
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Solution V depends on i?, which is denoted by v^. By compactness of we 
can find a subsequence converging to v such that 

^ ( 5 . 28 ) 

When Ce < S, then v € Ss, and {; is a solution of equation (5.5). 


5.7. Uniqueness. Because of the decay of v at infinity, we can obtain the unique¬ 
ness of V by the comparison principle as follows: 

Suppose that vi and V 2 are two solutions of (5.5). The difference w = vi — V 2 
satisfies the same equation and boundary conditions on the 2 ; 2 -axis, and w = 0 on 
the zi-axis. 

For any small positive constant r, we let R to be large enough such that |?c| < t 
on the cutoff boundary L^. Similar to (5.14), we employ Theorem A.l to obtain 
|w| < T in Q^. Let i? —?► oo and t oo, we conclude that w = 0, which implies 
the uniqueness. 


5.8. Determination of as a fixed point. We define a map Q : by 

Q{w) = V for any w = (p — 

where the closed set E^ is defined in (4.11). We employ the Schauder fixed point 
theorem to prove the existence of a fixed point for Q. That is, we need to verify 
the following facts: 

(i) E ,5 is a compact and convex set in a Banach space B; 

(ii) Q : Ei —>• Ea is continuous in B. 

Choose the Banach space B as defined in (4.12). Then E^ is compact and convex 
in B. 

For the continuity of Q, we make the following contradiction argument. Let 
wq,w'^ G Ei and ic " — >■ wq in B. Then u ” = Q(w^) in E^ and vq = Q{wq) in E^. 
We want to prove that v" —>■ vq in B. 

Assume that u" vq. Then there exist cq > 0 and a subsequence {i;"''} such 
that Iji;”'' — nolle > co- Since {n"'=} C E^ is compact in B, we can find another 
subsequence, again denoted by {n"''}, converging to some vi € E^. Then vg and 
vi satisfy the same equation (5.5), where ip = ip°° wg, which contradicts with the 
uniqueness of solutions for (5.5). Therefore, Q is continuous in B. 

Thus, we have a fixed point v for Q, which gives a solution ip = p°° + v for the 
nonlinear equation (4.2). The solution is unique by applying the same comparison 
principle as for the linear equation. 

Therefore, for given {A, U~ , b), we have determined tp. 


6 . Determination of the Entropy Function A as a Fixed Point. In this 
section, we employ the implicit function theorem to prove the existence of a fixed 
point A. 

6.1. Setup for the implicit function theorem for A. Through the shock polar, 
we can determine the values of U at O, and hence A(0) = At is fixed, depending 
on the values of U-{0) and b'{0). Then we solve (4.8) to obtain a unique solution 
A = h{z, (/?, V(/?) that defines the iteration map. To complete the proof, we need to 
prove that the iteration map exists and has a fixed point by the implicit function 
theorem. 

In order to employ the implicit function theorem, we need to set up a Banach 
space for A. To realize this, we perform the following normalization for (A, U~, b). 
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Let Aq = Define a smooth cutoff function x on [0, oo) such that 


X(s) = 


1, 0<s<l, 

0, s > 2. 


Let uj = U — Uq and ^ = b — Bq. Set 


A(0) 


( 6 . 1 ) 


where t is a function determined by the Rankine-Hugoniot conditions (3.6)-(3.9). 
Set X = A — wt with 

wt{z2) = - A+)xiz2). ( 6 . 2 ) 

Then A(0) = 0. 

Given (A,a;,/i), we can compute A — wt = J{A — wt) that defines the iteration 
map, by constructing a map X = A — wt =V{X,uj,p). 

We will prove that equation V{X, w, /r) — A = 0 is solvable for A, given parameters 
{ujyfj-) near (0,0). This is obtained by the implicit function theorem. 


6.2. Properties of the operator V. We first define some Banach spaces for 
operator V. Set 

X = {A : A(0) = 0,||A|U <oo} (6.3) 


with 


ll^llv = 


II \ ||(-a);{o} 
II^II1,q:;(1+/3);(0.oo)’ 


(6.4) 


Y = {uj : ||w||y < oo} (6.5) 

with 

l|w||v = ||w||2,a;(/3+l);n_ (6.6) 

for a vector-valued function w, and 

Z = {/r : /r(0) = 0, ||/r||z < 00 } (6.7) 

with 

ll/^b ^ ll/^lll,a;(/3);(0.oo)- (6-8) 

Clearly, X,Y, and Z are Banach spaces. Operator P is a map from X x Y x Z 
to X. 

We now define a linear operator D\'P{X,U!, n) and show that it is the partial 
differential of V with respect to A. When no confusion arises, we may drop the 
variables (A,w,/r) in D\V{X,uj, ^). 

We divide the proof into four steps. 


1 . Definition of a linear operator D\V{X,lu, fi). Given dX G X, we solve the 
following equation for Sip: 


Y, ( E 4iSip),^+btSip + dtSX)^^=0, 

(6.9) 

i=l,2 i^l,2 


with boundary conditions: 


Sip\z^^o = 0, 

(6.10) 

( E +c^<5(^)Ui=o = 0, 

2 ^1,2 

(6.11) 
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where 


= 5v:.., (2,^/5, V(/?), c^ = g^(z,(/3, V(/?). 

Once we have known (5 (/j, we define 


(5A = DxViX, w, /r)(5A) := ^ (6-12) 

z^l,2 


where 


V(p), = h^{z,ip,yip). 

It is easy to see that 5A(0) = 0. Then D\'P{X,uj, fi) is a linear operator from X 
to X. 

2. Show that D\V{X,u}, is the partial differential of V with respect to X at 
{X,u;,p). 

For fixed (w, p), we need to estimate 'P{X+6X, uj, p)—'P{X, w, p)—D\'P{X, lo, p){6X) 
to be o(i5A). 

For A, we define ip by following the definition of V, i.e., we solve the following 
equation, an alternative form from (4.2): 


{Xf{z,wt+X,p,Vp))^_=0, (6.13) 

i=1.2 

with boundary conditions (4.7) and (4.9). 

For A + AA, the corresponding potential (p satisfies 

^ (M*(z,u;t + A + AA,<p,V<^))^. =0, (6.14) 

i=l,2 


with the same boundary conditions (4.7) and (4.9). 

Taking the difference of equations (6.13) and (6.14) leads to the following equa¬ 
tion: 


H -^)+ ^A)^. = 0, (6.15) 

2 ^ 1,2 j^l,2 

with boundary conditions: 


{^-^)\z-2=0 = 0 , 

( X! ~ 1^1=0 = 

2 = 1,2 


(6.16) 

(6.17) 
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where 

= f {z,wt + X + s6X,(p + s{(p - ip),\7{ip + s{(p - ip)))ds, 

Jo 

6^ = [ ]Vf^{z,wt + X + s6X,(fi + s{(f - ip),\7{ip + s{(f - ip)))ds, 

Jo 

df^ = ( M\[z,wt + X +sSX,!^ + s{(p - ^p),S/{^p + s{(p - ip)))ds, 

Jo 

+ V5),V((^ + s((^-V3)))ds, 

Jo 

= / 9ip{'Z‘:<f + s{ip-if),V{if + s{ip-if)))ds. 

Jo 

Take the difference of (6.15) and (6.9), and \et v \= (p — — 5p. Then we have 

H H 

2^1,2 i^l,2 

= - E (E («'/ - 
2^1,2 3^1,2 

= (6.18) 
with boundary conditions: 

^1^2=0 = 0, 

^ uf u,. + ^ + (c^ - for zi = 0. 

2^1,2 2^1,2 

Since (p — ip = 0{dX) and 6ip = 0{6X), we conclude that and 2 i^i ~ 
i^i^)iS(p)zi + (c^ — c^^)S(p are o(i5A). Thus, we have 

V = o{SX). 


Then 


r{X + (5A, w, p) - V{X, to, p) - DxV{X, w, pi){5X) 

■=\-x-Tx 

= /i(z, <^, V(^) - /i(z, (j), V0) - ^ ei{Sp)z, - egSp 

2 ^1,2 

= E + E + (4^ - eo)'^‘P, 


where 


ef^ = / K., (z, (p + s((p- (p),V((p + s((p - ip))) ds, 
Jo 

/ h^{z,ip +s{(p-ip),V{ip +s{(p-ip)))ds. 
Jo 


pSX _ 

Cn — 


Therefore, we conclude that A — A — i5A = o((5A). Thus, D\'P{X,ijj, fi) is the partial 
differential of V with respect to A at (A, w, p,). 
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3. Continuity ofV and D\V. It suffices to show the continuity of V at any point 
{X*,u!*, fj,*) near (0,0,0) by a contradiction argument, since the same argument 
applies to DxV. 

Assume that there exists a sequence (A”,a;",/r") —>■ (A*,w*,/r*) in X x Y x Z, 
while IIA” — A*||x > cq > 0. Using the compactness of A” in || • H^' and the 
compactness of (/j" in || • ||i with a' < a, we find a subsequence {uk} such that 
Xrik converges to some A** in || • ||a/ and converges to some (p** in || • ||i,a'. Now 
we see that p* and p** satisfy the same equation (4.2) with the same boundary 
conditions. By the uniqueness of solutions for (4.2), we conclude that <p* = ip**. 
This implies that A* = A**. However, by assumption, ||A** — A*||x > cq > 0. This 
leads to a contradiction. Therefore, V is continuous. 

4. Show that, at the background state {X,uj,p) = (0,0,0), £>>,7^(0,0,0) — I is an 
isomorphism. 

When {X,u},p.) = (0,0,0), we solve for 6p: 

Y, {M\dX + = 0 in g, (6.19) 

2 = 1,2 

with the boundary conditions: 

E 1-1=0 = 0’ (6.20) 

2 = 1,2 

= 0 , ( 6 . 21 ) 

where {M\, pi) are the corresponding {Ma,M^^ evaluated at the back¬ 

ground state (C/1,1/0 ,(^q). 

Then we have 

DxViSX) ■.= sx=Y 

i=l,2 

where hi := h^^, evaluated at the background state. 

We rewrite the system in the following way: Let 


6p = Sp + m / SX(s)ds. 

Mf ^ ^ Jo 


Then (6.19)-(6.21) become 


^ Mj(d^),.,. =0 inQ, 

(6.22) 

Y = ™52^A, 

2 = 1,2 

(6.23) 

dp\z.,=o = 0. 

(6.24) 

Then 

(5A = y] hi{Sp)zi{0,Z2) - h2mdX. 

i=l,2 

Equations (6.23) and (6.25) give rise to 

(6.25) 

{h^i - gih2){.5p):zi = hSX. 

(6.26) 


TRANSONIC FLOWS PAST CURVED WEDGES 


27 


Noticing that i5A(0) = 0 and g 2 hi — 51/12 ^ 0, the boundary conditions (6.26) and 
(6.24) are compatible to guarantee the unique solution of (6.22) for arbitrary SX. 
This implies that I — DxV is onto. 

When 6X — D\'P{5X) = 0, (6.25) becomes 


h^{5lf)z^{0,Z2) = {I+ h2m)5X. (6.27) 

i=l,2 

Canceling 5X in (6.23) and (6.27) implies 

(51 + m{gih 2 - 52 ^ 1 )) ((5(/7)zi + g 2 {S^)z.z = 0. (6.28) 

The oblique boundary condition (6.28) above is nondegenerate, since 52 7 ^ 0. 
Therefore, solving equation (6.22) with boundary conditions (6.24) and (6.28) leads 
to dip = 0. Using (6.23), we conclude that 5X = 0, which implies that I — D\V is 
one-to-one. Thus, I — DxV is an isomorphism. 

Therefore, given U~ and 6, operator V has a fixed point A, which determines 
A = wt + X. 

With A from A, we obtain a unique potential ip. Using the Hodograph transfor¬ 
mation (4.1), we can solve for (j), with estimate 




2.a;(/3,0);Da 


< Ce. 


The subsonic flow U can be expressed by a smooth function of {X/(j),A,B), i.e. 
U(y) =F(Al(y 2 ), 5 ( 2 / 2 ), V<(.(y)). Also, = F [A, B .V cj)^). Then 

U-U°°= [ Fs^^{A,B,V(t)°° + s{V(l)-V(l)°°))ds{V(j)-V(l)°°) 

Jo 

= G{A,B,yv)yv. 


Based on the expression above, we can conclude 

Since the Lagrangian coordinate transformation is bi-Lipschitz, we can change the 
y-coordinates back to the x-coordinates and complete the existence part of Theorem 

2 . 1 . 


7. Decay of the Solution to the Asymptotic State in the Physical Coor¬ 
dinates. Now we determine the decay of the solution to the asymptotic state U°° 
in the x-coordinates. We divide the proof into four steps. 


1. For the fixed point established in §5.7, estimate (5.28) implies 

if-if^G Ecs. 

Then the change of variables from the z-coordinates to y-coordinates yields 


'^0|l2,a;(/3);R+ 




where is the subsonic region defined in (3.12). 

2. From (6.4), (7.1), and Step 4 in §6.2, we have 

Then, from §5.1, we obtain that, for f/°° = U°°{y 2 ) = {u°°, p°°){y 2 ), 


(7.1) 


(7.2) 


(7.3) 
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and 

u°° • (sin0o, — cos^o) = 0. (7.4) 

3. Since X 2 = ^(y), we now estimate (j){y) — tan^oj/i, which is X 2 — tanO^xi. 
From §5.1 and (7.1), 

(j)°° - tan6»o?/i = l{y 2 ), 

so that 

\\{(j) - tan6»o yi) - ^(l/ 2 )|l 27 J.[^“o)p^ < Ce. (7.5) 

In particular, this implies that, for each ?/2 > 0, 


l'{y 2 )= lim dy^{(j}{y)-ta.neoyi) = lim (l)y^{y). 

yi^oo yi^oo 

By (3.13), ^y, = By (7.2), 

^ I — Ce> — ^ , for small e > 0. 

Po ^10 2pq UiQ 

Therefore, there exists C > 0 such that, for any y 2 > 0, 

^ < I'M < c. 

Since i/’(0) = 0, we conclude that /(O) = wq, where wq is the limit of b 
§5.1). 

Furthermore, by (7.5)-(7.6), and (5.1) with (7.2), we have 

11/ _ 2/2 ||*,(-l-«);{0} „ 

Then there exists g : [wq, oo) —>■ [0, oo) with g{wo) = 0 such that g = l~^: 

9 {Ky 2 ))=y 2 on(0,oo) 


(7.6) 


bo (see 


(7.7) 


and 

^ < g'is) < C, 

so that g{-) satisfies (7.7). In order to composite g with (j) — tanPoyi, we need 
to extend g to a larger domain. Notice that (j)y^ > — ^ ^ implies that (f) is an 

2Pq 

increasing function in 2 / 2 - Therefore, condition (3.26) and the assumption (2.11) on 
b yields 

4 > — tan 00 yi > —£• 

Thus, we will extend I to domain [—wi, 00 ) so that g is defined on [—£, c»). By (7.7), 
we know that l'{y 2 ) > + — Ce > —-r—r, for small e. Take wi = and 

Po “10 ^Po “10 

we use the same extension map as in Lemma 6.37 in [17] to extend I to [—r(;i,oo) 
as follows. The extended function is still denoted by 1. Let 

h{s) = /(s) - 

PoKo 


For s G [—■u;i,0), define 

2 

h{s) = ^ Cih{-s/i), 
i=l 

where Ci, C 2 are constants determined by the algebraic equations 
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Then define 


l{s) = h{s) + 

Po ^10 


for s G [—^1,0). The extended function I satisfies 

ii7 y2 II* 


ny 2 ) > 


+ + lll,a;(/3);[-u;i,oo) 

Po “10 

1 


< Ce, 


2Po"'“io 

Let —W 2 = l{—wi) and we know 


for y 2 G [-wi,oo). 


W 2 =—l{0) + f l'{s)ds 

^ —Wi 


Wi 


>£ {|wo| < £ by (2.11)). 

Thus, g = is defined on [—£, oo). Therefore, by (7.5), we have 
\\9{<Piy) - tan6»oyi) - y2||l.a;(/3,0);Ds < Ce. 

In the same manner as for /, we can extend C/°° to domain [—wijOo) such that 
\\U°° — f^o"||o,c«;(l+,fl);[-™i,oo) < Cs. 

We define 

= U°°{g{s)), for s e [-£,oo). 

Then we employ (7.3) to obtain 

ll^°°(</>(y) - tan6»o?/i) - C/°“(2/2)||o,a;(/3,i);Ds <Ce. (7.8) 

4. Next, we use that the change of variables y —>■ x is globally bi-Lipschitz, 
which follows from (7.1) and (3.13) that implies the Jacobian: 


if £ is small, by (7.3). 

We also note that, in the y-coordinates, (7.1) implies for U°° = U°°{y 2 ) that 


\\u-u°-\\ 


{-a)-,aw 

l,a;(/3,l);Da 


< Ce. 


Then, changing the variables from y to x (which is bi-Lipschitz) and using (7.8), 
we obtain that, in the x-coordinates with X 2 = <(>(y). 


||17 - V°^{x 2 - tan6»o a;i)||o,a;(,3,i);as < Ce. 

This completes the proof for the decay of solution U{x) to the asymptotic state 
U°°. 
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8. Stability of Solutions. In this section, we prove that the subsonic solutions 
are stable under small perturbations of the incoming flows and the wedges as stated 
in Theorem 2.1. We modify operator V into V as follows: 

We first modify the definitions of the spaces in (6.4)-(6.7) in §6.2 by discarding 
the constraints: 


^ = {A : ||A||x < oo}, Z = {/i : ||/r||z < oo}, 

where norms || • ||x and || • \\z are defined in (6.4) and ( 6 . 8 ), respectively. We still 
use the same space Y and the related norm as in (6.5) and ( 6 . 6 ). 

Let w = U~ — Uq , X = A — ^(J", and p, = 6 — 69 . Given (A,u;,/r), we define 
P(A, uj, in the same way as for V by the end of §5, except that we do not restrict 
the value of ^(0) by (6.1). The restriction for A(0) is essential for the isomorphism 
of D\V. To prove the stability, we need to eliminate this restriction so that the 
differentiability in w can be achieved in a larger space. 

Equation (4.2) can be written as 

^ M*(C/_((^,Z 2 ),A 1 + + A,V<^),, = 0 . ( 8 . 1 ) 

i=l,2 


Given ((5A, 5uj, S^) G X xY x Z, define D'P{SX, Suj, 6^) in the following way: We 
solve the following equation for Sip: 


2^1,2 i^l,2 


with the boundary conditions: 

I _ i5/i(b) - S(I>-(J), 0) 

( 6 , 0 )’ 

( X! 1^1=0 = 0, 

2^1,2 


( 8 . 2 ) 

(8.3) 

(8.4) 


where 

“y- = +A,V?,V(/5), 

bi = mIj-{U~{(p,z 2),A+ +X,if,Vif) ■ {U~)y,{if,Z2), 
Z2),A^ + A,(^, Vv?), 
ft = ]^u-(U~{ip,Z2),A+ + A,(/?, V(/?), 

{U~{<p,Z2),V(p), 

= gU-iU~iyy,Z2),'^ip) ■ iU~)y^{ip,Z2), 

= 9u-(.U~{ip,z2),Vip). 


Then we define 

SX = DV{X,uj, p,){SX, Xuj, A/i) := '^2 ^Hb‘p)zi + e^Sip + • Slu, 

i=l,2 


(8.5) 
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where 

= K.SU-{ip, Z2},V<p), 

So = hu-(iy'~(‘P,Z2),V<p) ■ (U-)y,(<p,Z2), 

W'l = hu-{U~{(p,Z2),^if). 

Following the same estimates as in §6.2, we can verify that DV is the differential 
ofV. _ 

Let TZ := V — I. In §6.2, given (w, fi), we can find the fixed point for A. Therefore, 
A is a function of denoted by A(a;,/i). Therefore, we have 

7?.(A(a;,/i), w,/x) = 0. (8.6) 

Suppose that there is another parameter (w,/i) so that 

7^(A(d;,/x),a;,/2) = 0. (8.7) 

Taking the difference of equations (8.6) and (8.7), we have 

DxTl{X,uj,^i){X - A) + L>(^^_^)(d; - a;,/2 - /x) + o(A - A,d) - a;,/2 - ^) = 0, 

where A = A(w,/x). Since D\TZ{X,U!, fi) is an isomorphism near the background 
state, by inverting D\TZ^ we obtain 

A - A = - w,/x - ^) + o(A - A, w - w,/x - /x). 

Therefore, we obtain the following inequality: 

||A - A||x < ^(llw - u)\\y + ll/x - /x||z), (8.8) 

which implies the stability of the solutions depending on the perturbation of both 
the incoming flows and wedge boundaries. 

9. Remarks on the Transonic Shock Problem when Uq Is on Arc TH. 

In this case, uf > 0 and > 0 in the boundary condition (5.10), which makes a 
significant difference from the case when Uq is on arc TS. Such a difference may 
affect the estimates, hence the smoothness of the solutions, in general. 

In particular, one may not expect a solution for case TS is (7^’“; it is generically 
only in (7“. 

For example, in the first quadrant, let domain OAB be the quarter of the unit 
disc. The oblique direction u = (-1,-1). Let u = sin(|), where (r, 0) are the 

polar coordinates. In OAB, u satishes the Laplace equation: 

Au = 0, (9.1) 

and the boundary conditions: ?x = 0 on OA, and V?x • // = 0 on OB. However, u is 
Holder continuous only in (7 2 . 

Therefore, it requires a further understanding of global features of the problem, 
especially the global relation between the regularity near the origin and the decay 
of solutions at infinity, to ensure the existence of a smooth solution, more regular 
than the Holder continuity. A different approach may be required to handle this 


case. 
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Appendix A. Two Comparison Principles. In this appendix, we establish two 
comparison principles. 

Suppose that is a bounded, connected, and open set in M". Define a uniformly 
elliptic operator 

L= X] <93:, ( ^ -\-bi{x.)) in Q 

2^1,2 i^l,2 


in the following sense: 

for any x S D and ^ G M", 

id=1.2 


where A is a positive constant. Assume that aij,bi G C^(fl) fl C(0). 


Theorem A.l. Suppose that v,w G fl C(0) satisfy 

(i) Lv > 0 and Lw < 0 in D; 

(ii) w > 0 in Cl. 

Then — achieves its positive maximum on the boundary: 


sup 

n 



< sup 
an 



Proof. Let 


V 


V 

5 

w 


By calculation, we have 




2E 


h. 


Li=l,2 


E 7 -> T ^ Lm Lv 

BiVxi H- V = —. 

i=r ,2 ^ ^ 


(A.1) 


(A.2) 


By assumption, we know that — < 0 and — > 0. Therefore, by the weak maxi- 
mum principle. Theorem 8.1 in [17], we conclude (A.l). □ 


Theorem A.2. Suppose that v,w G <7^(17) fl C(D) satisfy 

(i) Lv > Lw and Lw < 0 in Ll; 

(ii) w > 0 in fl. 
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Then — achieves its positive maximum on the boundary or no greater than 1 in fl: 


/ v\ 

< max < sup 1 

fu^\ A 

sup - 

1 — )>i r 

n ^wJ 

t an 

\ w J ) 


(A.3) 


Proof. Equation (A.2) implies 


^ ^ {aijVx^)xj A ^ ^ Bi^Xi — 


Lw ^ , Lv — Lw 

- 1-E +- 

w w 


Lw 


{1-V). 


(A.4) 


Assume that V achieves the maximum value M > 1 at some interior point Xq G fl. 
Then, by continuity of V, there exists a ball Bm = Br(xo) C fi such that 


sup V = sup V = M > 1, 
Bm 

E > 1 in Bm- 


Therefore, 

— (1-E)>0 inSM, 
w 

and (A.4) implies 

^ ^ {^ij^xi')xj A ^ ^ BiYxi A 0 in Bm- 

i,j=l,2 i=l,2 


By the strong maximum principle. Theorem 8.19 in [17], we conclude 

V = M in Bm- 


(A.5) 


This implies that 

^ ( {s^ij^xi]xj A ^ ] BiV^xi =0 in Bm: 

i,j—1,2 2—1,2 

which contradicts (A.5). This completes the proof. □ 


Appendix B. The Shock Polar. We consider the uniform constant transonic 
flows with horizontal incoming supersonic flows. We now employ the Rankine- 
Hugoniot conditions (3.6)-(3.9) to derive a criterion for different arcs TS and TH 
on the shock polar. 

Assume that U~ and U are constant supersonic and subsonic states, respectively. 
The shock-front is a straight line: yi = sy 2 - Let k = ^ and k~ = = 0. Then 

the Rankine-Hugoniot conditions (3.6)-(3.9) give rise to 


[—] 
pup 

\ui A 


= —ks, 


P 


- = -pks, 

pup 


uik = [pjs, 
-1 


IP 


[2^1(1 + ^'^+ ( 7 -l)pJ 


= 0 . 


(B.l) 


(B.2) 

(B.3) 

(B.4) 
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From (B.3), s = Replacing s in (B.l) and (B.2), we obtain 


[- ][p] + uik'^ = 0 , 

[ui + ^^1 [p] + uipk'^ = 0. 

pwi 


(B.5) 

(B. 6 ) 


From equations (B.4)-(B.6), we can solve p,ui, and k in terms of p. Regarding 
(p, Ml, k) as functions of p, we differentiate (B.4)-(B.6) with respect to p to obtain 


BX = f, 


(B.7) 


where 


B = 


\p] 

fc2- J4, 

2miA: 

pUtT’ 

P“i ’ 

[P] 

-- 5 

^ p U\ ’ 


2puik 

IV 

-Mi(fc^ -h 1), 

—u\k 

(7-l)p2> 


X — {Pp: (Wl)p, kp) , 

pui pui pui 

We solve equation (B.7) for kp to obtain 

_ pCp 

" " Co ’ 


uik^ 


( 7 - 1)P 


where 


Cp = [p] (c^ + (7 - - 7 M 1 ) + (7 - l)pg^U 2 + 

Co = ulu2P^ ((7 + 1)P + (7 - 1)P") • 


LpMi 


(B. 8 ) 


p^C^Uiq^, (B.9) 
(B.IO) 


Notice that Cq > 0. Then, when state U belongs to TS, we find that kp > 0, which 
is equivalent to 

Cp < 0. (B.ll) 

On TH, we obtain that Cp > 0. 
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